VEHICLE CONTROL METHOD AND VEHICLE CONTROL APPARATUS 



Cross-Reference to Related Application 

This application claims priority under 35 USC 1 19 from 
Japanese Patent Application Nos. 2003-24177 and 2003-385973, the 
disclosures of which are incorporated by reference herein. 

BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a vehicle control method and a 
vehicle control apparatus. Specifically, the present invention relates to 
a method and an apparatus for controlling a vehicle which can control 
a tire force, that is, a force generated between a tire and a road surface, 
by independently controlling each steering angle and braking force of 
four wheels, and a vehicle which can control a tire force by 
independently controlling steering angles and braking forces of front 
wheels and those of rear wheels. 
Description of the Related Art 

A technique for separately controlling four wheels and steering 
angles thereof has been known, as disclosed in Japanese Patent 
Application Laid-Open (JP-A) No. 2001-322557. According to this 
related art, when a vehicle turns, the vehicle is steered so that a 
steering angle of each wheel becomes 90° with respect to the center of 
the turn of the vehicle. Further, when a defective wheel is found, 
steering and braking are controlled so that the friction force of the 
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defective wheel is decreased. This related art does not disclose 
cooperation of steering and braking, or cooperation of steering and 
driving. For example, regarding the steering angle, a fixed value is 
outputted as a target value regardless of braking or driving. 

However, in an actual vehicle, critical friction is generated 
between a wheel and a road surface, and lateral force may be decreased 
by increasing braking force. Therefore, cooperation between steering 
and braking or between steering and driving is indispensable in order to 
use the friction force between the wheel and the road surface as 
efficiently as possible. However, the above related art does not mention 
the cooperation at all, such as adjusting steering according to braking 
or driving. As a result, a problem arises in that the force generated by 
the wheel cannot be optimized. 

SUMMARY OF THE INVENTION 
An object of the present invention is to solve the 
aforementioned problem. 

In order to achieve the above object, a first aspect of the 
present invention may comprise: calculating a physical quantity which 
relates to a tire force of each wheel and optimizes an effective road 
friction of each wheel, based on a target resultant force to be applied to 
a vehicle body in order to obtain vehicle body motion that a driver 
desires, and a constraint including as parameters a magnitude of a 
critical friction circle of each wheel; calculating, based on the 
calculated physical quantity relating to the tire force of each wheel, a 
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first control variable for controlling at least one of braking force and 
driving force of each wheel, or a second control variable for controlling 
the first control variable and a steering angle of each wheel; and 
controlling (A) the at least one of braking force and driving force of 
each wheel based on the first control variable, or controlling (A) the at 
least one of braking force and driving force of each wheel and (B) the 
steering angle of each wheel, based on the first and second control 
variables. 

The constraint may be represented by a formula indicating 
that no resultant force is generated in a direction orthogonal to the 
direction of the target resultant force, and a formula indicating that a 
moment around the center of gravity of the vehicle is equal to a desired 
moment. Further, the constraint may be represented by formulae, the 
number of which is less than that of wheels, or a linearized formula. 

The target resultant force may be represented by a secondary 
performance function including a magnitude of the critical friction 
circle of each wheel and the physical quantity relating to the tire force 
of each wheel. In such a case, a physical quantity relating to the tire 
force of each wheel, which physical quantity satisfies a first 
approximation formula of the formula defining the constraint and 
optimizes the secondary performance function, can be calculated as 
the physical quantity which relates to the tire force of each wheel and 
optimizes the effective road friction of each wheel. 

Alternatively, a physical quantity relating to the tire force of 
each wheel, which physical quantity satisfies a first approximation 
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formula of the formula defining the constraint and optimizes the 
secondary performance function, is calculated as an initial value. The 
formula defining the constraint is linearized by using the calculated 
initial value. Then a physical quantity relating to the tire force of each 
wheel, which physical quantity satisfies the linearized formula of the 
constraint and optimizes the secondary performance function, is 
calculated as an approximate solution. The physical quantity which 
relates to the tire force of each wheel and optimizes the effective road 
friction of each wheel can be calculated by using the calculated 
approximate solution as the initial value to repeat the linearization of 
the formula defining the constraint and the calculation of the 
approximate solution. 

The formula defining the constraint may be linearized by 
Taylor expansion around the initial value or the approximate solution. 

The physical quantity relating to the tire force may be a 
direction of the tire force. The effective road friction of each wheel, the 
calculated direction of the tire force of each wheel, and the magnitude 
of the critical friction circle of each wheel can be used to calculate a 
slip angle based on a brush model, and the calculated slip angle can be 
used to calculate the second control variable based on a vehicle 
motion model. 

The magnitude of the critical friction circle of each wheel can 
be determined based on an estimate or a virtual value of a coefficient 
of friction p of each wheel and load of each wheel. 
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The direction of the tire force which optimizes the effective 
road friction of each wheel may be one of: a direction of the tire force 
which uniformly minimizes the effective road friction of each wheel; a 
direction of the tire force which makes the effective road friction of the 
front wheel differ from that of the rear wheel; and a direction of the tire 
force which makes the magnitude of the tire force of each wheel 
proportional to the load of the wheel. 

When the magnitude of the tire force proportional to the load 
of the wheel cannot be obtained because each wheel has a different \jl 
with respect to a road surface, the magnitude of the critical friction 
circle may be used as the magnitude of the tire force for a wheel having 
a small p, and the magnitude of the tire force which minimizes the 
effective road friction may be used for a wheel having a large p. 

The steering angle may be controlled so as to be the same for 
the right and left wheels. The effective road friction may be 
represented by the magnitude of the target resultant force relative to 
the magnitude of a critical resultant force obtained from the 
magnitude of the critical friction circle of each wheel. 

The direction of the tire force which is generated by each wheel 
may be defined as a value, that is the sum of products, which are 
calculated for all other wheels, of a distance from the position of an 
object wheel to the position of the other wheel in the direction of the 
resultant force, and the magnitude of the critical friction circle of the 
other wheel, with the direction of the resultant force acting on the 
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vehicle body as the resultant force of the tire forces of the respective 
wheels being used as a reference. 

A second aspect of the present invention may comprise: target 
resultant force calculating means for calculating a target resultant 
force to be applied to a vehicle body in order to obtain a vehicle body 
motion that a driver desires; critical friction circle estimating means 
for estimating the magnitude of a critical friction circle of each wheel; 
tire force calculating means for calculating a physical quantity which 
relates to a tire force of each wheel and optimizes an effective road 
friction of each wheel, based on the target resultant force and a 
constraint including as parameters the magnitude of the critical 
friction circle of each wheel; control variable calculating means for 
calculating, based on the calculated physical quantity relating to the 
tire force of each wheel, a first control variable for controlling at least 
one of braking force and driving force of each wheel, or a second 
control variable for controlling the first control variable and a steering 
angle of each wheel; and control means for controlling (A) the at least 
one of braking force and driving force of each wheel based on the first 
control variable, or controlling (A) the at least one of braking force and 
driving force of each wheel and (B) the steering angle of each wheel 
based on the first and second control variables. 

Further, a third aspect of the present invention may comprise: 
target resultant force calculating means for calculating a target 
resultant force to be applied to a vehicle body in order to obtain a 
vehicle body motion that a driver desires; critical friction circle 
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estimating means for estimating a magnitude of a critical friction 
circle of each wheel; critical resultant force estimating means for 
estimating a critical resultant force based on the magnitude of the 
critical friction circle of each wheel estimated by the critical friction 
circle estimating means; effective road friction setting means for 
setting a ratio of the target resultant force to the critical resultant 
force as an effective road friction; magnitude of tire force setting 
means for setting a magnitude of a tire force used at each wheel, which 
tire force is obtained by multiplying the magnitude of the critical 
friction circle of each wheel by the effective road friction; direction of 
tire force setting means for setting a direction of the tire force 
generated by each wheel based on a value, that is a sum of products, 
which are calculated for all other wheels, of a distance from the 
position of an object wheel to the position of the other wheel in a 
direction of the resultant force, and the magnitude of the critical 
friction circle of the other wheel, with the direction of the resultant 
force acting on the vehicle body as the resultant force generated by the 
tire force of each wheel being used as a reference; and control means 
for controlling a steering angle of each wheel and at least one of 
braking force and driving force of each wheel based on the magnitude 
and direction of the tire force which have been set. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a schematic view of a vehicle motion model. 
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Fig. 2 is a schematic view of a coordinate system 
corresponding to a resultant force generated in the vehicle motion 
model in Fig. 1. 

Fig. 3 is a diagram of a critical friction circle. 

Fig. 4 is a block diagram of a first specific structure according 
to an embodiment of the present invention. 

Fig. 5 is a diagram showing the results of simulation, 
indicating a comparison of vehicle accelerations when, in the first 
specific structure, the brake is applied to a vehicle which advances 
straight such that one side of the vehicle travels on a snowy road 
surface and the other side of the vehicle travels on a dry road surface. 

Fig. 6A is a diagram showing calculated directions of tire 
forces of respective wheels, and the number of times calculations are 
repeated, in the case of front-wheel and rear-wheel steering. Fig. 6B is 
a diagram showing calculated directions of tire forces of respective 
wheels, and the number of times calculations are repeated, in the case 
of four-wheel independent steering. 

Figs. 7A and 7B are diagrams showing vectors of the tire forces 
of the respective wheels and the steering angles of the respective 
wheels after convergence. 

Fig. 8 is a diagram showing tire force characteristics. 

Fig. 9 is a diagram showing the results of simulation, 
indicating a comparison of critical vehicle accelerations when the 
brake is applied to a vehicle which is turning such that one side of the 



8 



vehicle travels on a wet road surface and the other side of the vehicle 
travels on a dry road surface. 

Fig. 1 OA is a diagram showing calculated directions of tire 
forces of respective wheels, and the number of times calculations are 
repeated, in the case of front-wheel and rear-wheel steering. Fig. 10B 
is a diagram showing calculated directions of tire forces of respective 
wheels, and the number of times calculations are repeated, in the case 
of four-wheel independent steering. 

Figs. 1 1A and 1 IB are diagrams showing vectors of the tire 
forces of the respective wheels and the steering angles of the respective 
wheels after convergence. 

Fig. 12 is a diagram showing comparison between tire 
characteristics derived based on a brush model, which theoretically 
describes the tire force, and actual tire characteristics. 

Fig. 13 is a diagram showing vectors of the tire forces of 
respective wheels and the steering angles of the respective wheels after 
convergence, when the brake is applied to a vehicle which is turning 
on a road surface of uniform \i. 

Fig. 14 is a block diagram showing a second specific structure 
according to the embodiment of the present invention. 



DETAILED DESCRIPTION OF THE INVENTION 
A preferred embodiment of the present invention will be 
described in detail below with reference to the drawings. First, a 
principle of cooperative control of steering and braking in a vehicle 
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which enables independent steering and braking of four wheels, and a 
principle of cooperative control of steering and driving in a vehicle 
which enables independent steering and driving of four wheels will be 
described. 

In a motion model of a four-wheel-drive vehicle shown in Fig. 
1, a direction 9 in which a force is applied to a vehicle body as a 
resultant force, which is the sum of forces generated by four wheels, in 
order to obtain a vehicle body motion that a driver desires, and the 
magnitude (radius) F { of a critical friction circle of each wheel (wherein 
i = 1 to 4) are known, (i = 1 represents a left front wheel, i = 2 
represents a right front wheel, i = 3 represents a left rear wheel, and i 
= 4 represents a right rear wheel.) In this case, the direction of the tire 
force which maximizes the resultant force, namely, acceleration (or 
deceleration) generated by the vehicle body, while securing a desired 
yaw moment, is determined for each wheel. The direction of the tire 
force of each wheel is represented by an angle q i? which is formed by 
the direction of the resultant force and the direction of a force 
generated by a single wheel (i.e., the tire force of each wheel). 

The critical friction circle is a circle representing the limit in 
which motion performance of the vehicle can be controlled without 
losing grip of a tire. The magnitude of the critical friction circle 
represents the maximum friction force of a tire generated between the 
wheel and a road surface, and is determined based on an estimate or a 
virtual value of p (friction coefficient) and load of each wheel. The 
friction force of the tire is formed by a force in a direction in which the 
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vehicle advances (driving force or braking force) and a friction force in 
a lateral direction (rightward or leftward). The friction force in either 
direction becomes 100% of the magnitude of the critical friction circle, 
or corresponds thereto, the friction force in the other direction 
becomes zero. The braking force is applied in a direction opposite to 
the direction in which the driving force is generated. The range of the 
friction force can be represented as a substantial circular shape when 
represented as a vector diagram as shown in Fig. 3. For this reason, 
the range of the friction force is referred to as the critical friction 
circle. 

When the motion model of the four-wheel-drive vehicle in Fig. 
1 is converted into a coordinate shown in Fig. 2 with an X-axis being a 
direction of the resultant force and a Y-axis being a direction 
perpendicular to the X-axis, the position of each tire, i.e., (x, y) = (b i? 
a { ), can be represented by the following formulae (1) to (8). 
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a x = 
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— - cos 0 — L f sin 0 
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(1) 


a 2 = 
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— — - cos 0 — L f sin 0 
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a 3 = 
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— L cos 0+ L r sin 0 
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a 4 = 
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— — L cos 0+ L r sin 0 
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b 2 = 
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— — ^ sin 0+ L f cos 0 
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(6) 


b 3 = 


T 

— - sin 0— L r cos0 
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(7) 


b 4 = 


T 

— — 11 sin 0— L r cos0 
2 


(8) 



In the formulae, T f is an interval between front wheels, Tr is an 
interval between rear wheels, L f is a distance between the center of 
gravity of the vehicle and a middle point of the interval between the 
front wheels, and L r is a distance between the center of gravity of the 
vehicle and a middle point of the interval between the rear wheels, ^ 
and bj represent a distance from the X-axis and the Y-axis, 
respectively. 

Further, assuming that M z represents a yaw moment (desired 
moment) to be generated around the center of gravity of the vehicle at 
this time, constraints represented by the following formulae (9) and 
(10) exist for the angle q { , which represents the direction of the tire 
force of each wheel. 
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F : sin q 2 + F 2 sin q 2 + F 3 sin q 3 + F 4 sin q 4 = 0 ( 9 ) 



— a l F 1 cos q x — a 2 F 2 cos q 2 — a 3 F 3 cos q 3 — a 4 F 4 cos q 4 

+ b^sin qj+ b 2 F 2 sin q 2 + b 3 F 3 sin q 3 + b 4 F 4 sin q 4 =M J , 

(10) 

The formula (9) represents a constraint indicating that 
resultant force is not generated in the Y-axis direction, namely, a 
direction orthogonal to the direction of the resultant force. The 
formula (10) represents a constraint indicating that the moment 
around the center of gravity of the vehicle is a desired yaw moment M z . 
The number of the formulae representing these constraints is less 
than the number of the wheels. 

Therefore, the problem of maximizing the resultant force, 
namely, the problem of maximizing an effective road friction, becomes 
a problem of determining the angle qi which satisfies the constraints of 
the formulae (9) and (10) and maximizes the sum (resultant force) J of 
forces in the X-axis direction represented by the following formula 
(11). 

J = Fj cos q x + F 2 cos q 2 + F 3 cos Q 3 + F 4 cos q 4 (11) 

This problem can be solved as a nonlinear optimization 
problem by numerical calculation of convergence, which will be 
described later. Alternatively, a solution can be derived from 
approximation described as follows. 
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First, the formulae (9) and (10) representing the constraints 
are subjected to primary approximation to obtain the following 
formulae (12) and (13). 

Fj Q x + F 2 q 2 + F 3 q 3 + F 4 q 4 = 0 (12) 

q x + b 2 F 2 q 2 + b 3 F 3 b 3 + b 4 F 4 q 4 

= M Z + a 1 F 1 + a 2 F 2 + a 3 F 3 + a 4 F 4 (13) 

The formula (11) is further subjected to secondary 
approximation to obtain a secondary performance function 
represented by the following formula (14). 

J = F -^i 2+ p 2 -^ 2+ F 3-^ q 3 2 + F 4-^U 2 (14) 

Further, the problem of determining the angle q { which 
maximizes the secondary performance function represented by the 
formula (14) can be replaced with a problem of determining the angle 
q i which minimizes the sum (K/ 2) of negative terms in the formula (14), 
namely, a problem of determining the angle q { which minimizes a 
secondary performance function represented by a formula (15). 

Namely, the problem of maximizing the effective road friction 
becomes a problem of determining the angle q { which maximizes the 
secondary performance function that is represented by the formula 
(14) and includes the magnitude of the critical friction circle of each 
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wheel and the direction of the tire force of each wheel, or a problem of 
determining the angle q { which minimizes the secondary performance 
function that is represented by the following formula (15) and includes 
the magnitude of the critical friction circle of each wheel and the 
direction of the tire force of each wheel, namely, a problem of 
determining the angle q { which optimizes the effective road friction. 

K = F lQl 2 + F 2 q 2 2 + F 3 q 3 2 + F 4 q 4 2 (15) 

Variables in the formula (15) are changed by using the 
following formula (16). 

Pi^VFTq, (16) 

The formula (15) is represented by the following formula (17). 

K= Pl 2 +p 2 2 +p 3 2 +p 4 2 (17) 

As a result, the above problem turns out to be replaced with a 
problem of determining a value p { of the minimum Euclidean norm 
which satisfies the following formulae (18) and (19). 
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^ Pi +Vf7p2+#7p3+#7p4= 0 (18) 

b 1A /F7Pi+ b 2 VF7p 2 + b 3 7F7p 3 + b 4 VF7p 4 
= M Z + ajF^ a 2 F 2 + a 3 F 3 + a 4 F 4 (19) 



Pi can be solved as in the following formula (20). 
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F, VF 3 VF 4 
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M z + a^H- a 2 F 2 + a 3 F 3 + a 4 F 4 
"c. 



M 2 + a,F,+ a 2 F 2 + a 3 F 3 + a 4 F 4 
bjC,F,+ b 2 c 2 F 2 + b 3 c 3 F 3 + b 4 c 4 F^ 



L C 4 



(2 0) 



The formula (20) satisfies the following formulae (21) to (24). 



c x = (b-b 2 ) F 2 + (b-b 3 ) F 3 + (bf-b^ F 4 (2 1) 

c 2 = (b 2 -b x ) F x + (b 2 -b 3 ) F 3 + (b 2 -b 4 ) F 4 (2 2) 

c 3 = (b 3 -b x ) Fj+ (b 3 -b 2 ) F 2 + (b 3 -b 4 ) F 4 (2 3) 

c 4 = (b 4 -b,) F t + (b 4 -b 2 ) F 2 + (b 4 -b 3 ) F 3 (2 4) 
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Further, "diag" in the formula (20) represents a diagonal 
matrix with values in the parentheses being diagonal elements. 
represents a pseudoinverse matrix. When A is a long full rank matrix, 
the pseudoinverse matrix of A can be calculated by the following 
formula (25). 

A + =(A- A T )"'A T (2 5) 

As will be described later, the angle q { derived from the above 
formulae (20) to (24) is directly used as the direction of the tire force of 
each wheel when the effective road friction y is 1, in integrated control 
of the steering angle and the braking force or the steering angle and 
the driving force. 

Further, the formulae (21) to (24) respectively represent the 
sum of products, which are calculated for all other wheels, of the 
distance from the position of the wheel where the angle q { is 
determined to the position of the other wheel in the direction of the 
resultant force of the vehicle body (e.g., (bj-bs), (b!-b 3 ), (bj-b^ for a 
wheel with i = 1), and the magnitude of the critical friction circle of the 
other wheel. Thus, these formulae represent that the angle q { formed 
by the direction of the tire force of each wheel and the direction of the 
resultant force of the vehicle body is proportional to the sum of the 
products, which are calculated for all other wheels, of the distance 
from the position of the wheel where the angle q { is determined to the 
position of the other wheel in the direction of the resultant force of the 
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vehicle body, and the magnitude of the critical friction circle of the 
other wheel. 

Furthermore, a { or b 4 representing the position of each wheel is 
a function of the direction 9 of the resultant force of the vehicle body. 
Therefore, the angle q 4 formed by the direction of the tire force of each 
wheel and the direction of the resultant force of the vehicle body can 
be represented as a function of the direction of the resultant force of 
the vehicle body and the magnitude of the critical friction force of each 
wheel. 

The angle q { derived as described above can also be used as an 
initial value for the formulae (9) to (11) for convergence calculation of 
the nonlinear optimization. Generally, in the nonlinear optimization 
problem, convergence in the convergence calculation tends to 
accelerate by using a value near the optimal value as an initial value. 
For this reason, solutions of the formulae (20) to (24) can be used as 
the initial values in the nonlinear optimization problem to increase the 
efficiency of calculation. 

In the nonlinear optimization problem mentioned herein, 
Taylor expansion around approximate solutions derived from the 
formulae (20) to (24) for the formulae (9) and (10) is carried out to 
derive solutions therefrom. Further, Taylor expansion around the 
solutions for the formulae (9) and (10) is repeated to derive solutions. 
Consequently, approximate solutions of high accuracy are derived. 

First, the angle qi is derived, as an initial value q i0 , from the 
first Taylor expansion, using the formulae (20) to (24), for the formulae 
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(9) and (10). The following formulae (26) and (27) are derived from 
Taylor expansion around the initial value q i0 . 



F 1 {sinqj 0 +(Qi _ Qi o )cosq! 0 H" F 2 (sinQ 2 o+fe - ^2 o )cosq 2 0 } 

+ F 3 {sinq 3 0 +(q 3 - Q 3 0 )cosQ 3 0 }+ F 4 {sinq 4 0 +(q 4 - ^4 o )cosq 4 0 }= 0 

(2 6) 

- a 1 F 1 {cosq! 0 _ (Qi _ Qi o )sinQi 0 }~ a 2 F 2 {cosq 2 0 -(q 2 _ Q 2 o ) sin Q 2 o } 
- a 3 F 3 {cosq 3 o _ (q 3 _ ^3 o )sinq 3 0 }- a 4 F 4 {cosq 4 o - ^ - OU o )sinq 4 0 } 

+ b 1 F 1 [smQ 1 0 +(Qi- Qi 0 )cosQ 1 0 }+ b 2 F 2 {sinq 2 0 +(q 2 - q 2 0 )cosq 2 0 } 
+ b 3 F 3 {sinq 30 +(q 3 - q 3 o )cosq 3 0 }+ b 4 F 4 {sinq 4 0 +(q 4 ~ q 4 o)cosq 40 } 

= M z 

(2 7) 



A solution which satisfies these formulae (26) and (27) and 
also minimizes the formula (15) is derived from the following formula 
(28), in which the pseudoinverse matrix described above is used. 
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1 1 1 1 

V^i" _ 

TFTCaiSinqj 0 + bj cosq! 0 ) ^/F7(a 2 sinq 2 0 + b 2 cosq 20 ) 

VF7cosq 30 V^4" cosC l4o 
^F7(a 3 sinq 3 0 + b 3 cosq 3 0 ) jF^fa 4 sin Q4 o + b 4 cosq 4 0 ) 

4 " 

L F i(Q> o cosQj Q-sinqj 0 ) 
M z + IF, {(a,+ b,q, 0 )cosq, 0 +(a,q, 0 ~ b, )sinq, 0 } 

i — 1 

(2 8) 



Q 2 
Q 3 



= diag 



Next, solutions which satisfy the formula obtained by the 
Taylor expansion around approximate solutions of the formula (28) 
and minimize the formula (15) are derived by using the pseudoinverse 
matrix, as described above. Subsequently, Taylor expansion and 
derivation of approximate solutions are repeated a predetermined 
number of times to derive an angle q { of high accuracy. 

The accuracy of optimization can be improved by repeatedly 
performing calculations with the formula (28) being used as a 
recurrence formula, namely, by using the q i? which has been 
calculated in the previous step, as q i0 in the next step and repeating 
Taylor expansion and the calculation of q^ 

Next, cooperative control of steering and braking and of 
steering and driving before a limit will be described. Cooperation of a 
steering sj'stem and a braking system, and cooperation of a steering 
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system and a driving system for improving critical motion performance 
(i.e., the force generated by the vehicle body) have been described 
above. A cooperation method for extending the control rules obtained 
above to a range before the limit and maximizing grip of each wheel will 
be described below. 

When the effective road friction of each wheel is expressed as 
Y, constraints of a force generated by the vehicle body in a lateral 
direction and a yaw direction are represented by the following 
formulae, which are similar to the formulae (9) and (10), respectively. 

yF x sinq 1 +YF 2 s i n Q2 + YF 3 sinq 3 +YF 4 sinq 4 = 0 

— a^Fjeosq!— a 2 yF 2 cosq 2 — a 3 YF 3 cosq 3 — a 4 yF 4 cosq 4 

+ b 1 yF 1 sinq x + b 2 yF 2 sinq 2 + b 3 YF 3 sinq 3 + b 4 yF 4 sinq 4 =M z 

Namely, the constraints are represented by the following 
formulae (29) and (30). 

F x sinq^ F 2 sinq 2 + F 3 sinq 3 + F 4 sinq 4 = 0 (29) 

— a^jCosQi— a 2 F 2 cosq 2 — a 3 F 3 cosq 3 — a 4 F 4 cosq 4 

M 

+ b 1 F 1 sinq^ b 9 F 2 sinq 9 + b 3 F 3 sinq 3 + b 4 F 4 sinq 4 = — - (30) 

Y 

Further, the magnitude of the force generated by the vehicle 
body is represented by the following formula as a constraint. 
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yFj cosq^yFg cosq 2 +Y F 3 cosq 3 +vF 4 cosq 4 = F 

Namely, the constraint is represented by the following formula 

(31). 

p 

F x cosq^ F 2 cosq 2 + F 3 cosq 3 + F 4 cosq 4 =-j- (31) 

Therefore, the cooperation method for uniformly maximizing 
the grip of each wheel becomes a problem of determining the angle ^ 
which satisfies the formulae (29) to (31) and minimizes the effective 
road friction y. Further, when F * 0, the problem can be considered as 
a problem of determining the angle q { which satisfies the above 
formula (29) and the following formula (32) obtained by organizing the 
formulae (30) and (31), and maximizes the following formula (33). 



(3 2) 
(3 3) 

When an approximate solution is derived in the same way as 
described above, the formulae (29) and (32) are represented by the 



b 1 F 1 sinq^ b 2 F 2 sinq 2 + b 3 F 3 sinq 3 + b 4 F 4 sinq 4 
a x +— - Fjcosq^ 
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following formulae (34) and (35), respectively, after primary 
approximation. 



F 2 q 2 + F 3 q 3 + F 4 q 4 = 0 



(3 4) 



b 1 F 1 q 1 + b 2 F 2 q 2 + b 3 F 3 q 3 + b 4 F 4 q 4 = 



a x + — * 



M ^ 



f 



F 2 + 



M. ^ 



3 F 



F 3 + a 4 + 



M ^ 



(3 5) 



Further, the formula (33) corresponds with the secondary 
performance function of the formula (14) by secondary approximation. 
For this reason, as described above, variables are changed by using 
the formula (16), and the problem is replaced with a problem of 
determining a value p ; of the minimum Euclidean norm which satisfies 
the following formulae (36) and (37). 



/f7p 1 +7f7p 2 +Vf7p 3 + a /f7p 4 = o 



(3 6) 



b 1 #7Pi+ b 2V F 7P 2 + b 3 7F7p 3 + b 4 VF~P 4 = 



' . M ^ 



Fi+ a 2 + 



M. > 



a Q + 



M ^ 



F 3 + 



a A +— i 



(3 7) 



The value p ; can be solved by the following formula (38). 
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q 4 . 



=diag 



=diag 



V^s" _ 

111 



Pi 
p 2 
p 3 

l_p 4 J 



V^l V"^3 



Vf7 Vf7 #7 #: 



_bi#7 b 2 VF7 b a#7 b 4#7_ 

0 



(F.+ F 2 + F 3 + F 4 >f a.F^ a 2 F 2 + a 3 F 3 + a 4 F, 



_ F 



(F,+ F 2 + F 3 + F 4 >f ai F x + a 2 F 2 + a 3 F 3 + a 4 F, 



b^F^ b 2 c 2 F 2 + b 3 c 3 F 3 + b 4 c 4 F, 



(3 8) 



The formula (38) satisfies the following formulae (39) to (42). 



c 1 =(b-b 2 )F 2 +(b 1 -b3)F3+(b 1 -b 4 )F 4 (3 9) 

c 2 =(b 2 -b 1 )F 1 +(b 2 -b3)F 3 +(b 2 -b 4 )F 4 (4 0) 

c 3 =(b 3 - b, jF.+fe- b 2 )F 2 +(b 3 - b 4 )F 4 (41) 

c 4 =(b 4 - b x )Fi+(b 4 - b 2 )F 2 +(b 4 - b 3 )F 3 (42) 



As described above, "diag" represents a diagonal matrix, and 
represents a pseudoinverse matrix. The angle qi derived from the 
formulae (38) to (42) may be directly used in integrated control of 
steering and braking or steering and driving, as the direction of the 
tire force of each wheel, or may be used as an initial value for the 



formulae (29) to (31) for convergence calculation of the nonlinear 
optimization. In the same way as described above, a recurrence 
formula for determining an approximate solution of high accuracy by 
repeated calculation based on the Taylor expansion is given as the 
following formula (43). 



: diag 



yf^2 ^f^4 _ 



^fF~ 2 ~ cosq 2 o 



jjF^fa sind 1 0 + b : cosq! 0 ) V^T( a 2 sinC l2 o + b 2 cosQ 2 0 ) 

7F7cosq 30 7f7 cosq 4o 
-jF^(a 3 sinq 3 0 + b 3 cosq 3 0 ) (a 4 sinq 4 0 + b 4 cosq 4 0 ) 

4 

LFi^i ocosqj Q-sinq; 0 ) 



4 

I 

i = 1 



IF, j a s + b,q, o + cosqi 0 +(a,q, 0 - b.jsinq, 0 



(4 3) 



The effective road friction y is calculated based on the formula 
(44) in which the angle qi thus derived is used. Namely, the effective 
road friction y is represented as a ratio of a target resultant force (i.e., 
force of the vehicle body) F to a critical resultant force obtained from 
the magnitude of the critical friction circle of each wheel (i.e., 
resultant force in the direction of the target resultant force, which 
resultant force is the sum of critical forces obtained from the 
magnitude of the critical friction circle of each wheel). 
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Fj cosq^ F 2 cosq 2 + F 3 cosq 3 + F 4 cosq. 



(4 4) 



Moreover, braking or driving force of each wheel is derived 
from the following formula (45) by using the effective road friction y, 
and the magnitude F i of the critical friction circle of each wheel, and 
the direction (qj + 0) of the tire force of each wheel. yF { represents the 
magnitude of the tire force. 

F x r^FiCOS^+e) (4 5) 

Further, lateral force applied to each wheel is derived from the 
following formula (46). 

F y ^YF^infe+e) (4 6) 

The steering angle of each wheel is calculated based on, for 
example, a brush model and a vehicle motion model. The brush model 
is a model which describes a characteristic of the tire force based on a 
theoretical formula. Assuming that the tire force is generated in 
accordance with the brush model, a slip angle pi can be determined by 
the following formula (47), using the magnitude of the critical friction 
circle of each wheel, the effective road friction y, and the direction of 
the tire force of each wheel {q { + 0). 
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P,= tan 



_/ K s -k,sin(q,+ e) 



(4 7) 



K„ 1 - kjCosCqs+e) 



The slip angle $ { satisfies the following formula (48). 




1 -(1 - Y )§ 



J 



(4 8) 



In the above formulae (47) and (48), K s represents driving 
stiffness, and K p represents cornering stiffness. 



the slip angle based on the vehicle motion model. Namely, the steering 
angle can be calculated as in the following formulae (49) to (52) by 
using yaw angle speed r 0 and a vehicle body slip angle (3 0 , which are 
calculated as target vehicle motion state variables from vehicle speed 
v, the steering angle, accelerator stroke, braking effort, and the like, 
to (3 4 represent the slip angles of the respective wheels in the formula 



Further, a steering angle 8 ; of each wheel is calculated from 



(47). 



5 2 = / ? 0 + ~ r O _ ^2 
V 

_ L f 



(4 9) 



(5 0) 



(5 1) 



(5 2) 
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When cooperation of steering control and braking and driving 
control is carried out, braking force and driving force are controlled 
based on first control variables, which are the braking force and 
driving force determined by the above formula (45). The steering 
angle, namely, the direction of the tire force, is controlled based on 
second control variables, which are the steering angles determined in 
the formulae (49) to (52). Alternatively, either the braking force or the 
driving force may be controlled while the direction of the tire force is 
controlled. 

When the cooperation of the steering control and the braking 
and driving control is carried out based on this type of control, the 
effective road friction y of each wheel can be uniformly minimized all 
the time, and motion performance allowing the greatest robust against 
disturbance such as a road surface or cross wind can be obtained. 

Further, when the resultant force to be applied to the vehicle 
body is maximized, the braking force and the driving force of each 
wheel can be determined by the formula (45) when the effective road 
friction y of each wheel is 1, and the steering angles of the respective 
wheels can be determined by the formulae (49) to (52) when the 
effective road friction y of each wheel in the formula (48) is 1. 

The braking force, the driving force and the steering angle of 
each wheel determined above are used as the control variables to 
control the cooperation of the driving force and the steering angle of 
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the vehicle or the cooperation of the braking force and the steering 
angle. 

Only the braking force and the driving force may be calculated 
to control the driving force, the braking force, or both the braking force 
and the driving force. Namely, only the magnitude of the tire force may 
be controlled without controlling the steering angle. 

Next, a case will be described in which the cooperative control 
is applied to a normal four-wheel vehicle and carried out by using the 
same steering angle for the right and left wheels. In the case of the 
conventional four-wheel vehicle having the same steering angle for the 
right and left wheels, constraints represented by the following 
formulae (53) and (54) and indicating that the right and left wheels 
have the same slip angle are further added. 



tan 



- 1 



K 1 - k^cos^+e) 



tan" 



^K Sj -k 3 sin(q 3 +e) 
K p 1 — k 3 cos(q 3 +0) 



= tan 1 



K s . ~ k 2 sin(q 2 +6) > 
K p 1 — k 2 cos(q 2 +0) 

K s m - k 4 sin(q 4 +6) 
K p 1 - k 4 cos(q 4 +e) 



(5 3) 
(5 4) 



These constraints can be organized into the following formulae 
(55) and (56). 



k 2 sin(q 2 +6)- kiSin^+e)- k : k 2 sin(q,- Q 1 )= 0 
k 4 sin(q 4 + 0) — k 3 sin(q 3 + 9) — k 3 k 4 sin(q 4 — q 3 ) = 0 



(5 5) 
(5 6) 



After primary approximation, these constraints are 
represented by the following formulae (57) and (58). 



k^cosS— k 2 )q!+ k 2 (cos0- k 1 )q 2 =(k 1 — k 2 )sin6 
k 3 (cos6- k 4 )q 3 + k 4 (cos0- k 3 )q 4 =(k 3 - k 4 )sin0 



(5 7) 
(5 8) 



As described above, constraints such as the following 
formulae (59) and (60) are used when the effective road friction y is 1, 
namely, when the force of the vehicle body is maximized. 

F : sinq x + F 2 sinq 2 + F 3 sinq 3 + F 4 sinq 4 = 0 (59) 

— a^jcosq^ ^2^2 cosq 2 ~ a 3F 3 cosq 3 — a 4 F 4 cosq 4 

+ b l F l sinq^ b 2 F 2 sinq 2 + b 3 F 3 sinq 3 + b 4 F 4 sinq 4 =M 7 (60) 

Therefore, an approximate solution is derived from the 
following formula (61). 
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k 2 (cos6— k x ) 
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— k 3 (cos6 



b 4 F 4 
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k 4 (cos6- k 3 ) 



- 1 
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M z + a.F^ a 2 F 2 + a 3 F 3 + a 4 F 4 
(k : - k 2 )sin0 
(k 3 — k 4 )sin0 



(6 1) 
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Further, when the force F of the vehicle body is applied with 
the effective road friction y being 1 or less, the formula (32) is used 
instead of the formula (60). Furthermore, the effective road friction y 
is determined by the following formula. 

= F 

Y F 1 +F 2 +F 3 +F 4 

An approximate solution is derived from the following formula 
(62) by using the effective road friction y determined above. 







Fi 


F 2 F 3 




F 4 


Q 2 




bxF, 


b 2 F 2 b 3 F 3 
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+ ai F,+ a 2 F 2 + a 3 F 3 + a 4 F 4 




(6 2) 






(k- 


k 2 )sin£? 










(k- 


k 2 )sin# 







The braking force, the driving force and the steering angle at 
this time are calculated in the formulae (45) and (47) to (52). However, 
the same value is obtained in the formulae (46) to (49) for the steering 
angles of the right and left wheels. The solutions of the formulae (61) 
and (62) are derived by primary approximation. As described above, 
these solutions may be used as the initial values for the corresponding 
nonlinear equation for numerical calculation, and the control based 
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on a solution thus obtained may be carried out. When the repeated 
calculations using the above-described Taylor expansion are applied 
to the formula (61), an approximate solution is represented by the 
following formula (63). 











q 3 




q 4 . 





F, cosq, o F 2 cosq 2() 
a,F, sinQ, 0 + b,F, cosq, 0 a 2 F 2 sinq 20 + b 2 F 2 cosq 20 

- k, {cos(q, 0 + d) ~ k 2 cos(q 2 0 - q, 0 )} k 2 {cos(q 2 0 + 6) - k, cos(q 2 „- Q, 0 )} 
0 0 

F 3 cosq 3 o F 4 cosq 4 0 

a 3 F 3 sinq 30 + b 3 F 3 cosq 3 o a 4 F 4 sinq 4 0 + b 4 F 4 cosq 4 0 

0 0 

k 3 {cos(q 3 0 + d) - k 4 cos(q 4 0 - q 3 0 )} k 4 {cos(q 4 0 + 8) - k 3 cos(q 4 „- q 3 0 )}. 

4 

52 F sCQi ocosq, Q-sinq; 0 ) 



M 7 .+ £F i {(a i + b,q, 0 )cosq, 0 +(a,Q, 0 _ b, )sinqi 0 } 



i = 1 



(6 3) 



The formula (63) satisfies the following formula. 



dj= kj {sin^ 0 + d)~ q, 0 cos^ 0 + d)) 

- k 2 {sin(q 20 +(9)- q 2 0 cos(q 2 0 + d)} 

+ kjkojsinCqa 0 - q x o)-(q2 o~ Qi o) cos (q 2 o _ Qi o)J 
d 2 = k 3 {sin(q 3O +0)- q 3 0 cos(q 3 0 + #)} 

- k 4 {sin(q 4 0 + d)~Q 40 cos(q 4 0 + e)} 

+ k 3 k 4 {sin(q 40 - Q 30 )-(q. 4 o- Q 3 o)cos(q 4 0 ~ q 30 )} 



q j0 : value obtained in the previous step 
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Further, when the repeated calculations using the Taylor 
expansion are applied to the formula (62), an approximate solution is 
represented by the following formula (64). 
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(6 4) 



Next, the first specific structure of the present embodiment 
using the above principle will be described based on Fig. 4. As shown 
in the drawing, the present embodiment includes a target resultant 
force calculating unit 18 for calculating the magnitude and direction 
of a target resultant force; a critical friction circle estimating unit 20 
for estimating the magnitude (radius) of a critical friction circle of each 
wheel; a direction of tire force and effective road friction calculating 
unit 25 for calculating the direction of force generated by each wheel 
and the effective road friction based on the magnitude and direction of 
the target resultant force and the magnitude of the critical friction 
force of each wheel; a tire force setting unit 28 for setting the 
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magnitude of the tire force generated by each wheel; and a control unit 
30 connected to a cooperative braking and driving apparatus. 

The target resultant force calculating unit 18 calculates, from 
the steering angle, the vehicle speed, the accelerator stroke, the 
braking effort, and the like, the magnitude and direction of a resultant 
force and a yaw moment to be applied to the vehicle body in order to 
obtain a vehicle body motion that the driver desires. The magnitude 
and direction of the resultant force and the yaw moment can be 
determined by, for example, calculating the magnitude and direction 
of the resultant force and the yaw moment which are necessary to 
approximate to zero deviations from a yaw angle speed, which is a 
target vehicle motion state variable set in accordance with the 
operation of the driver, and from a measured value (or an estimate) of 
the slip angle of the vehicle body. 

The critical friction circle estimating unit 20 estimates the 
magnitude of the critical friction circle of each wheel based on SAT 
(self-aligning torque) or a wheel speed. 

The direction of the tire force and effective road friction 
calculating unit 25 calculates the direction of the tire force and the 
effective road friction of each wheel based on the above constraints 
including as parameters the magnitude and direction of the target 
resultant force and the magnitude of the critical friction circle of each 
wheel. In this case, the direction of the tire force may be calculated so 
that the effective road friction of each wheel is uniformly minimized, or 
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may be calculated so that the magnitude of the tire force of each wheel 
is minimized in accordance with a load distribution ratio. 

When the magnitude of the tire force proportional to a wheel 
load cannot be calculated because each wheel has a different \i with 
respect to the road surface, the magnitude of the critical friction circle 
may be used as the magnitude of the tire force of the wheel having a 
small p, and the magnitude of the tire force obtained with the smallest 
effective road friction may be used for the wheel having a large p. 

As described above, the direction of the tire force of each wheel 
can be determined based on linear algebraic equations, the number of 
which is less than the number of wheels. Further, when the effective 
road friction of each wheel is uniformly minimized, as described above, 
the direction of the tire force can be determined based on the sum of 
products, which are calculated for all other wheels, of the distance 
from the position of an object wheel, the direction of the force of which 
is to be determined, to other wheel, which distance is in a direction in 
which a resultant force is applied to the vehicle body as the resultant 
force of the tire forces, and the magnitude of the critical friction circle 
of the other wheel. 

The tire force setting unit 28 calculates the steering angle, the 
braking force and the driving force of each wheel based on the 
direction of the tire force of each wheel for optimizing the effective road 
friction. 

The control unit 30 controls a steering actuator, a braking 
actuator and a driving actuator so that the steering angle, the braking 
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force, and the driving force calculated by the tire force setting unit 28 
are obtained. 

Fig. 5 is a diagram showing the results of calculation, by 
simulation, of the comparison of vehicle accelerations [G] when, in the 
first specific structure described above, the brake is applied to the 
vehicle which advances straight such that one side of the vehicle 
travels on a snowy road surface having a p of 0.3 and the other side of 
the vehicle travels on a dry road surface having a p of 1 .0. In this case, 
control based on the formulae (28) and (63) is carried out in four-wheel 
independent steering and front-wheel and rear-wheel steering, 
respectively, and it can be understood that braking force is increased 
by 8% in the four-wheel independent steering. 

Figs. 6A and 6B are diagrams showing the number of 
calculations repeated of the direction of the tire force, which 
calculations are based on the formulae (28) and (63), respectively. 
Figs. 7 A and 7B are diagrams each showing a vector of the tire force 
and a steering angle of each wheel after convergence. 

It can be understood from Figs. 6A and 6B that the direction of 
the force generated by each wheel converges in a second or third 
calculation. Further, it can be understood from Fig. 7B that, in the 
case of the four-wheel independent steering, a lateral force canceling a 
moment is generated at a wheel having a high p, which wheel has 
relatively large friction allowance. 

The front-wheel and rear-wheel steering shown in Fig. 7A 
adopts a control rule which uses an edge of the friction circle (k a in a 



36 



tire characteristics shown in Fig. 8). However, considering the use of 
a slip area (k b in the tire characteristics shown in Fig. 8) which is 
further than the edge, the direction of the tire force can be changed 
without changing the steering angle. 

In this case, the direction of the tire force of each wheel, which 
is the solution for the four-wheel independent steering, can be realized 
by a front-wheel and rear-wheel steering vehicle. The acceleration 
generated by the vehicle body which is equivalent to that of the four- 
wheel independent steering can be obtained by the following logic. 

First, a slip angle and a steering angle are determined based 
on the calculated value of the direction of the tire force of a wheel 
having a high ya (a control rule of the four-wheel independent steering). 
The wheel of high yi outputs a braking force and a driving force which 
are the same as those in the four-wheel independent steering. 

Next, a slip ratio for making the slip angles of the wheels of 
high \i correspond to each other and obtaining the direction of the tire 
force calculated from the control rule of the four-wheel independent 
steering is determined as follows. When a low yi and a high yi are 
represented by subscripts i and j, respectively, a condition in which 
the slip angles of the right and left wheels correspond to each other is 
represented by the following formula: 



d& s thar k i sin(q i +9) 
K 'l - kicostqi+e) 



= tan 



- 1 



^K s m ~ kj sin(q j + 8) 
^K p 1 -kjcosfej+e^ 



(6 5) 
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' K. 



In the above formula (65), q { and qj are solutions (directions of 
the tire force) derived by assuming the four-wheel independent 
steering. These solutions are used as follows to solve kj. 



3F J sra(q J +e) 
1 K s sin(q i +0)+ 3 F.sin(q - qj 



The slip ratio (in a longitudinal direction of the vehicle) is 
calculated as follows. 



k x .= k, cosfa+6)^ : \' ' . / (67) 

K s sin(q i +0 > )+ 3 F^in^— qj 



Further, the braking force or the driving force is given by the 
following formula. 

Fxi^FiCosCqi+O) (6 8) 

Fig. 9 is a diagram showing the results of calculating, by 
simulation, the comparison of critical accelerations [G] (the directions 
of the resultant force) during such turning and braking that the 
resultant force is generated in a direction of 9 = 120° on a wet road 
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surface having ap of 0.8 for outer wheels and a dry road surface 
having ap of 1.0 for inner wheels. 

"FRONT-WHEEL AND REAR-WHEEL STEERING" indicates an 
acceleration of the vehicle when the steering angles of the front and 
rear wheels and the braking forces of the respective four wheels are 
determined based on the formula (63) in the vehicle structured such 
that the right wheel and the left wheel are steered at the same angle. 
Further, "FOUR-WHEEL INDEPENDENT STEERING" indicates an 
acceleration of a vehicle when the direction of the force generated by- 
each wheel is determined based on the formula (28). It can be seen 
that the acceleration generated by the vehicle body under these 
conditions is increased by 5% in the four-wheel independent steering. 

Figs. 10A and 10B are diagrams showing the number of 
calculations repeated of the direction of the tire force of each wheel, 
which calculations are based on the formulae (28) and (63), 
respectively. Figs. 1 1A and 1 IB are diagrams each showing a vector of 
the tire force and an steering angle of each wheel after convergence. It 
can be seen that, in the case of the four-wheel independent steering, 
the steering angles of the outer wheels having a small p are smaller 
than those of the inner wheels. Cornering power (CP), dependency of 
lateral force on the load, a change in susgeometry, and the like are not 
considered in this simulation. 

In this simulation as well, the front-wheel and rear-wheel 
steering shown in Fig. 10A adopts the control rule which uses the edge 
of the friction circle (k a in the tire characteristic in Fig. 8). However, 
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considering the use of the slip area further than the edge (k b in the tire 
characteristics in Fig. 8), the direction of the tire force can be changed 
without changing the steering angle. Therefore, as in the case of the 
vehicle advancing straight such that the one side and the other side of 
the vehicle travel on different types of road surface, the front-wheel 
and rear-wheel steering vehicle, in which the right and left wheels are 
controlled so as to have the same steering angle, can also obtain a 
vehicle acceleration which is equivalent to that of the four-wheel-drive 
vehicle. As described above, in a critical traveling state (wherein the 
effective road friction = 1), the front-wheel and rear-wheel steering 
vehicle in which the right and left wheels have the same steering angle 
can also obtain the vehicle acceleration which is equivalent to that of 
the four-wheel independent steering by using the slip area further 
than the edge of the friction circle. This is the result of derivation 
based on the brush model which theoretically describes the tire force. 
However, in actual tire characteristics, tire force may decrease in the 
slip area as shown in Fig. 12. When the solution for the four-wheel 
independent steering is realized in the front-wheel and rear-wheel 
steering vehicle, the vehicle acceleration may decrease by a 
corresponding amount. 

Fig. 13 is a diagram showing solutions (vectors of forces 
generated by the respective wheels and the vehicle body, and the 
steering angles thereof) for the four-wheel-drive vehicle and the 
front-wheel and rear-wheel steering vehicle when the brake is applied 
to the vehicles during turning on a flat road surface having a p of 0.85. 
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There is no difference between the four-wheel independent steering 
and the front-wheel and rear-wheel steering on the road surface 
having a uniform ]a, and both the four-wheel independent steering and 
the front-wheel and rear-wheel steering obtain the solution of q { = 0. 

Considering that almost no difference is found between the 
four-wheel independent steering and the front-wheel and rear-wheel 
steering vehicle on the road surface having a uniform p, and that the 
direction q { of the force generated by each wheel is relatively small, the 
following control rule for the front-wheel and rear-wheel steering may 
be used in critical and pre-critical traveling states. 

First, in the pre-critical traveling state, the maximum value of 
the road surface \x of each wheel is used, and control when the road 
surface ]i is maximum and uniform is carried out. When the wheels do 
not have a uniform \i, the actual effective road friction does not 
become uniform, but the magnitude of the tire force of each wheel 
becomes proportional to the load distribution ratio, and the direction 
q { of the force generated by each wheel becomes relatively small. 
Therefore, well-balanced cooperation of steering and braking, and of 
steering and driving can be expected. 

The control rule for the four-wheel independent steering is 
used when the wheel having a low |i has reached its limit, namely, in 
the case of the following formula (69). 

F i Rea.<YF i (6 9) 
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Y is calculated from the solution qi of the recurrence formula 
(64) and calculated as follows. 

Y = (70) 

Fj cosqiH- F 2 cosq 2 + F 3 cosq 3 + F 4 cosq 4 

F i is the magnitude of the tire force when |i of the road surface 
is assumed to be large, and F iReal is an actual magnitude of the tire 
force of a wheel having a small p. yF { in the process of deriving the 
control rule for the four-wheel independent steering can be obtained 
as F iReal . For example, when wheels 1 and 3 (the left front wheel and 
the left rear wheel) have reached their limits, the constraints for the 
wheels are represented as follows. 

FiReai sinq 1 +YF 2 sinq 2 + F 3 R e a , sinq 3 +YF 4 sinq 4 = 0 (71) 

~ a 1 F 1Real cosq 1 - a 2Y F 2 cosq 2 

- a 3 F 3Rea i^osq 3 - a 4Y F 4 cosq 4 
+ b i F iRea isinq^ b 2 yF 2 sinq 2 

+ b 3 F 3Rea isinq 3 + b 4 vF 4 sinq 4 =M z (7 2) 

F 1Rea iCosq 1 +vF 2 cosq 2 + F 3Rea l cosq 3 +yF 4 cosq 4 = F (7 3) 

Thus, this becomes a problem of determining q { which satisfies 
the formulae (71) to (73) and minimizes y- While the nonlinear 
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optimization may be used, the following recurrence formula may also 
be used: 



<*4 



=diag 



e a 1 



e a 1 



VYo F 4_ 



V F lReal COS Q, o VY0 F 2 COSq 2 0 

V F iRea i( a i sinQi 0 + b x cosq! o ) VYo F 2 ( a 2 sinq 2 0 + b 2 c°sQ 2 0 ) 



3 R e a 1 
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V F 3Rea I C °SQ 3 0 VYo F 4 COSq 4 0 

(a 3 sinq 3 0 + b 3 cosq 3 0 ) VVo F 4 ( a 4 sinq 4 0 + b 4 cos Q4 o )_ 



(7 4) 



provided that 

d 3 =L F i Rea i^i ocosqi o-sinqj 0 )+ LYo^fe o ^sQ, o-sinQi 0 ) 

i = l,3 i = 2,4 

d 4 =M z +£F i Re a ,{(a,+ b,q, oJcosq, o+faq, 0 - b, )sinq i 0 } 



i=l,3 



+ I Yo F i K a i+ tVli o )cosq, 0 +(a,q, 0 _ b i )sinq, 0 } 

i=2,4 



and 



y=- 



i=l,3 



(7 5) 



/=2,4 



The subscript "0" in Yo anc * F i0 indicates a value calculated in 
the previous step. Based on q 4 thus derived, the braking force, the 
driving force and the steering angle can be derived from the formulae 
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(45) to (52) for the wheel of large p, and, for the wheel of small p, the 
braking force, the driving force and the steering angle can be derived 
from the formulae (65) to (68), considering that the slip of the wheel of 
large ]i is as follows. 



As described above, when the control rule for the vehicle 
motion according to the present embodiment is derived with the 
direction of the tire force of each wheel being used as the control 
variable, a nonlinear tire model does not need to be included in the 
optimal calculation, and the number of control variables is smaller, as 
compared to a conventional method in which the slip angle and the 
slip ratio of each wheel are used as the control variables. Therefore, 
the amount of calculation of the control rule is smaller than that of the 
conventional method, and the use of the control rule can be expanded 
to the four-wheel independent steering, which is a system of a higher 
degree of freedom. 

In the four-wheel-drive vehicle, the case has been described in 
which the control rule which equalizes and minimizes the effective 
road friction of each wheel is derived. This means maximizing friction 
allowance of each wheel in combinations of control of each wheel for 
obtaining desired force and moment. Therefore, improvement in 
safety and resistance to failure can be expected. 




(7 6) 
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Moreover, the case has been described in which the control 
rule which equalizes the effective road friction of each wheel and the 
control rule which makes the tire force of each wheel proportional to 
the load distribution are derived in the front-wheel and rear-wheel 
steering vehicle in which the right and left wheels have the same 
steering angle. These control rules correspond to each other when the 
front-wheel and rear-wheel steering vehicle travels on the road surface 
of uniform p, and substantially correspond to the solution of the 
above-described control rule for the four-wheel-drive steering. Thus, 
this means that, when the front-wheel and rear-wheel steering vehicle 
travels on the road surface of uniform p, these control rules maximize 
friction allowance of each wheel in combinations of control of each 
wheel for obtaining desired force and moment. Therefore, 
improvement in safety and resistance to failure can be expected. 

Further, in the control rule for making the tire force 
proportional to the load distribution in the front-wheel and rear-wheel 
steering vehicle in which the right and left wheels have the same 
steering angle, when the wheel of small p has reached its limit of 
friction, the direction of the force generated by each wheel, which 
direction is derived based on the control rule for the four-wheel-drive 
vehicle, is realized by the front-wheel and rear-wheel steering vehicle. 

Next, a second specific structure of the present embodiment 
using the above principle will be described with reference to Fig. 14. 

In the present embodiment, the present invention is applied to 
a vehicle having an electric power steering device mounted therein. As 
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shown in the drawing, the present embodiment is formed by a group of 
sensors 10 mounted in the vehicle and including a steering angle 
sensor for detecting a steering angle from a rotation angle of a steering 
shaft of the electric power steering device, an assist torque sensor for 
detecting power assist torque from electric current passing through 
the electric power steering device, a torque sensor for detecting 
steering torque, a vehicle speed sensor for detecting the speed of the 
vehicle, an accelerator stroke sensor for detecting an accelerator 
stroke, a braking effort sensor for detecting braking effort, and a wheel 
speed sensor for detecting the rotation speed of a wheel; a controller 
12 formed by microcomputers; and a cooperative braking and driving 
device 14 which is connected to the controller 12 and cooperatively 
controls the steering angle, and braking and driving. 

The controller 12 formed by the microcomputers is controlled 
so as to perform a plurality of functions in accordance with programs 
stored in advance. When the controller 12 is shown in a block 
functional diagram, the controller 12 is formed by a SAT estimating 
unit 16 for estimating the SAT; a target resultant force calculating unit 
18 for calculating the magnitude and direction of the target resultant 
force; a critical friction circle estimating unit 20 for estimating the 
magnitude (radius) of the critical friction circle of each wheel; a critical 
resultant force estimating unit 22 for estimating the magnitude of 
critical resultant force; an effective road friction setting unit 24; a 
magnitude of tire force setting unit 26 for setting the magnitude of tire 
force; a direction of tire force setting unit 28 for setting the direction of 
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the tire force of each wheel; and a control unit 30 connected to the 
cooperative braking and driving device 14. 

The SAT estimating unit 16 estimates SAT based on steering 
torque detected by the torque sensor and assist torque detected by the 
assist torque sensor. 

The target resultant force calculating unit 18 calculates the 
magnitude and direction of the target resultant force and the yaw 
moment M z to be applied to the vehicle body in order to obtain a 
vehicle body motion that the driver desires, from the steering angle, 
the vehicle speed, the accelerator stroke, the braking effort, and the 
like which have been detected by the respective sensors. 

For example, a resultant force and a yaw moment that are 
necessary to approximate to zero deviations from a yaw angle speed, 
which is a target vehicle motion state variable, and from a measured 
value or an estimate of the slip angle of the vehicle body, can be used 
as the magnitude and direction of the target resultant force and the 
yaw moment M 2 . 

The critical friction circle estimating unit 20 estimates the 
magnitude of the critical friction circle of each wheel based on the SAT 
estimated by the SAT estimating unit 16 or the wheel speed detected 
by the wheel speed sensor. 

The critical resultant force estimating unit 22 estimates the 
magnitude of the critical resultant force from the magnitude of the 
critical friction circle of each wheel estimated by the critical friction 
circle estimating unit 20. The estimation may be carried out by simply 
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regarding the sum of the critical friction forces of the respective wheels 
as the critical resultant force, or by regarding as the critical resultant 
force a value obtained by multiplying the sum of the critical friction 
forces of the respective wheels by a constant. Alternatively, the 
estimation may be carried out by first determining the angle q { based 
on the formulae (21) to (24) from a direction 0 of the target resultant 
force and a target yaw moment, and then calculating a critical 
resultant force J based on the formula (11). 

The effective road friction setting unit 24 sets as the effective 
road friction y the ratio of the magnitude of the target resultant force 
to that of the critical resultant force. However, when the target 
resultant force exceeds the critical resultant force, the effective road 
friction is set as 1. 

The magnitude of tire force setting unit 26 sets the magnitude 
of the tire force used at each wheel by multiplying the critical friction 
circle by the effective road friction. 

The direction of the tire force setting unit 28 uses the direction 
0 of the target resultant force and the magnitude yF { of the tire force of 
each wheel used, to determine the angle q; based on the following 
formula obtained by replacing F { in the formulae (20) to (24) by yF { , and 
outputs the magnitude yF { and the direction (q { + 0) of the tire force of 
each wheel: 



48 



Q 4 



= diag 



a/^I "7^2 V^3 _ 



#7 



^f 2 #7 vf; 
b 2 VF7 b 3 VF7 b 4A /F7. 

0 



M 



^+ (a 1 F 1 +a 2 F 2 +a 3 F 3 +a 4 F 4 ) 



M 



_ Y 



^+ ( ai F x + a 2 F 2 + a 3 F 3 + a 4 F 4 ) 



b^F^ b 2 c 2 F 2 + b 3 c 3 F 3 + b 4 c 4 F, 



(7 7) 



provided that 

c x = (b x -b 2 ) F 2 + (b 1 -b 3 ) F 3 + (b 2 -b 4 ) F 4 (7 8) 

c 2 = (b 2 -b x ) F x + (b 2 -b 3 ) F 3 + (b 2 -b 4 ) F 4 (7 9) 

c 3 = (b 3 -b x ) Fj+ (b 3 -b 2 ) F 2 + (b 3 -b 4 ) F 4 (8 0) 

c 4 = (b 4 -b x ) F x + (b 4 -b 2 ) F 2 + (b 4 -b 3 ) F 3 (8 1) 



The control unit 30 determines the steering angle and the 
braking force, or the steering angle and the driving force of each wheel 
based on the magnitude yF i and the direction (q { + 0) of the tire force of 
each wheel, and controls the steering device and the braking actuator, 
or the steering device and the driving actuator. The braking force or 
the driving force of each wheel can be derived as follows from the 
magnitude yF { and the direction (q { + 0) of the tire force of each wheel. 
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Fx ,= yF|Cos (qj+G) 



(8 2) 



Similarly, lateral force of each wheel can be derived from the 
following formula. 

F yi = Y F i sin(q i + 0) (8 3) 

The steering angle of each wheel can be calculated, for 
example, based on the brush model and the vehicle motion model. The 
brush model is the model which describes the characteristic of the tire 
force based on the theoretical formula. Assuming that the tire force is 
generated in accordance with the brush model, the slip angle $ { can be 
calculated by the following formula, using the critical friction circle F { , 
the effective road friction y, and the direction of the tire force (q d + 0): 



P^tan" 



J K s -kjSinCqi+B) 



K 1 - k iC os (q { + 0) _ 



(8 4) 



provided that 

k^iA [i- (i_ r )i/3 ] (8 5) 



In the above formulae, K s represents driving stiffness, and K p 
represents cornering stiffness. Further, the steering angle of each 
wheel is calculated from the slip angle fi { based on the vehicle motion 
model. Namely, the steering angle can be calculated as follows by 
using the yaw angle speed r 0 and the vehicle body slip angle p 0 , which 
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are calculated as the target vehicle motion state variables from the 
vehicle speed v, the steering angle, the accelerator stroke, the braking 
effort, and the like. 



5i = Po+^r 0 - 




(8 6) 




K 


(8 7) 


5 3 = Po-^i- 0 - 


\ 


(8 8) 


5 4 =(3 0 -^r 0 - 


K 


(8 9) 



When the cooperative control of steering and braking and of 
steering and driving is carried out based on this type of control, the 
effective road friction of each wheel can be made uniform all the time, 
and motion performance allowing the greatest robust against 
disturbance such as a road surface or cross wind can be exhibited. 

Further, in the present embodiment, the effective road friction 
can be independently set for the front wheel and the rear wheel. In 
this case, the following formula is constructed, wherein the effective 
road frictions of the front wheel and the rear wheel are \ f and y r 
respectively: 



Q 4 



_ M z + a x Yf F t + a 2 y f F 2 + a 3 Yr F 3 + a 4 Yr F 4 
b i c i Yf Fi+ b 2 c 2 y f F 2 + b 3 c 3 Yr F 3 + b 4 c 4 v r F 4 



C 2 



(10 0) 
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provided that 



c 1 =(b 1 -b 2 ) Yf F 2 +(b 1 -b 3 )Y r F3+(b 1 -b 4 )Y r F 4 (10 1) 

c 2 =(b 2 -b 1 ) Yf F 1 +(b 2 -b 3 )y r F 3 +(b 2 -b 4 )Y r F 4 (10 2) 

C3 == (b3-b 1 )Y f F 1 +(b3-b 2 )Y f F 2 +(b 3 -b 4 )Y r F 4 (10 3) 

c 4 = (b 4 - b, ) Yf F,+ (b 4 - b 2 ) Yf F 2 + (b 4 - b 3 ) Yr F 3 (10 4) 



As described above, the effective road friction is independently 
set for the front wheel and the rear wheel. Namely, the effective road 
frictions of the front wheel and the rear wheel are made to differ from 
each other. As a result, for example, by setting the effective road 
friction of the rear wheel smaller than that of the front wheel, friction 
allowance of the rear wheel is increased, thereby realizing motion 
control which is highly effective in suppressing spin and emphasizes 
safety. 

Further, the formulae (101) to (104) respectively represent the 
sum of products of the distance from the position of the object wheel to 
the position of the other wheel in the direction of the resultant force of 
the vehicle body, and the magnitude of the tire force to be outputted, 
which tire force is obtained in consideration of the effective road 
friction. As a result, the formulae (101) to (104) represent that the 
angle formed by the direction of the tire force of each wheel and the 
direction of the resultant force of the vehicle body is proportional to 
the sum of the products of the distance from the position of the object 
wheel to the position of the other wheel in the direction of the resultant 
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force of the vehicle body, and the magnitude of the tire force to be 
outputted. 

Further, considering that a f or bj is a function of the direction 
0 of the resultant force of the vehicle body, the formulae (101) to (104) 
represent that the angle formed by the direction of the force generated 
by each wheel and the direction of the resultant force of the vehicle 
body can be described as a function of the direction of the resultant 
force of the vehicle body and the magnitude of the tire force of each 
wheel to be outputted. 

Next, an embodiment will be described in which the present 
invention is applied to the normal four-wheel vehicle, which is 
controlled such that the steering angle is the same for the right and 
left wheels. 

In the case of the conventional four-wheel vehicle having the 
same steering angle for the right and left wheels, constraints 
represented by the following formulae (105) and (106) and indicating 
that the right and left wheels have the same slip angle are further 
added. These constraints can further be organized into the following 
formulae (107) and (108). 
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tan 



- 1 



K g . -lqsinfa+e) 1 
K p 1 - k lC os (q x + 6) 



= tan 1 



r K s ^ -k 2 sin(q 2 +e) "I 
K p 1 — k 2 cos (q 2 + 6 ) 



(10 5) 



tan 



- 1 



K Sa -k 3 sin(q 3 +e) 
1 -k 3 cos(q 3 +0) j 

^K Sj -k 4 sin(q 4 +e) ^ 
^K p 1 - k 4 cos (q 4 + 6 ) ; 



tan 



(10 6) 

k 2 sin (q 2 + 0 ) - k x sin (q x + 0 ) - k x k 2 sin (q 2 - Q 1 ) = 0 (10 7) 



k 4 sin (q 4 + 0 ) — k 3 sin (q 3 + 0 ) — k 3 k 4 sin (q 4 — q 3 )= 0 (10 8) 

After primary approximation, these constraints are 
represented by the following formulae (109) and (110). 



kj (cos 0 — k 2 )q x + k 2 (cos 0 — k x )q 2 

= (k x — k 2 )sin 0 

■k 3 (cos 0 — k 4 )q 3 + k 4 (cos 0 — k 3 )q 4 

= (k 3 — k 4 )sin 0 



(10 9) 



(110) 



Furthermore, these constraints are represented by using p { as 



follows. 
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k k 
-7= (cos 8 - k 2 )Pi+-t= (cos 8 - k x )p 2 

= (kj— k 2 )sin 6 



(111) 



k k 
3 (cos 6 - k 4 )p 3 +- ? = (cos 6 - k 3 )p. 



= (k 3 — k 4 )sin 0 



(112) 



By considering other constraints represented by the following 
formulae (113) and (114) when the vehicle travels with an effective 
road friction y, the angle formed by the direction of the resultant 
force and the force generated by a single wheel can be uniquely 
determined as shown in the following formula (115). 

VF7 Pi +Vf;p 2 +Vf;P3+#7p4 = 0 (113) 

b!#7Pi + b 2 #7p 2 + h 3 jF^v 3 + b 4A /F7p 4 

(^F^ a 2 F 2 + a 3 F 3 + a 4 F 4 ) (114) 

Y 
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Q 4 j 



= diag 



#7 

^^(cos0- k 2 ) 

#7 



0 



r(cos6— kj ) 



0 



_ r 3 



0 

r(cos6— k 4 ) 



M 



0 



^+a 1 F 1 +a 2 F 2 +a 3 F 3 + a 4 F 4 



(k x — k 2 )sin9 
(k 3 — k 4 )sin9 



- 1 



0 

r(cOS0— k 3 ) 



(115) 



The braking force, the driving force and the steering angle of 
each wheel at this time are calculated by the formulae (82) and (84) to 
(89). However, the same value is calculated in the formulae (87) and 
(89) for the steering angles of the right and left wheels. The value ^ 
obtained by formula (115) is calculated by primary approximation and 
may also be used as an initial value to numerically solve a 
corresponding nonlinear equation, and the control is carried out 
based on the solution of the equation. In this case, control of higher 
accuracy can be achieved. 

When the cooperation of the steering control and the braking 
and driving control is carried out based on this type of control, the 
effective road friction of each wheel can be made uniform all the time, 
and motion performance allowing the greatest robust against 
disturbance such as a road surface or cross wind can be exhibited. 
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Moreover, even when the same steering angle is used for the 
right and left wheels, the effective road friction of each wheel can be 
set independently. In this case, the angle q ; formed by the direction of 
the resultant force and the force generated by a single wheel is 
represented by the following formula. 



q 4 



: diag 



V F 1 Yf V F 2 Y f V F 3 Yr V F 4 Yr 



b iV F i Yf 



b 2 V F 2 Yf 



V F 3 Yr 



1 (cosg-k 2 ) 2 (cosg-k, ) 
V F > Yf V F 2 Yf 



0 



V F 4 Yr 



b W F 4 Yr 

0 



0 



k k 
- . 3 (cosfl-k 4 ) 4 (cosfl-k 3 ) 

V F 3 Yr V F 4 Yr 



M z +(ai F i K/+ a 2 F, K/ + a 3 F 3 >/ r + a 4 F 4 Yt ) 
(kj— k 2 ) sin0 
(k 3 — k 4 )sin0 



F > Yf 



'zYf 



b,F, 



F 3 Kr F 4 Kr 

b 3 F 3 Kr b 4 F 4 K/ 

k,(cos0-k 2 ) k 2 (cose-kj) 0 0 

0 0 — k 3 (cos0-k 4 ) k 4 (cos(9- k 3 )j 

0 

M^a.F, Yf + a 2 F 2 K/ + a 3 F 3 Kr + a 4 F 4 Kr ) 
(k,— k 2 )sin0 
(k 3 — k 4 )sin0 



(116) 



The present embodiment can also be applied to a vehicle 
structured such that only the front wheels or the rear wheels can be 
independently steered. For example, in a case of a vehicle in which 
each of the rear wheels can be independently steered, the angle q ; 
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formed by the direction of the resultant force and the force generated 
by a single wheel can be represented by the following formula. 



Qi 
Q 4 



= diag 



.V F lYf #777 V F 3Yr V F 4Yr_ 



V F lYf 



V F 2Yf 



# 3 Yr V F 4Yr 



^#777 



t> 2 V F 2Y f 



-— ^--(cosO- k 2 ) Jf 2 (cose- k x ) 



0 

M z +( a i F i Y f + a 2 F 2 Yf + a 3 F 3 Yr + a 4 F 4 Yr ) 
(k x — k 2 )sin0 



b 3 V F 3Y r 
0 



b 4 V F 4Yr 

0 

(117) 



The magnitude of the tire force in the present embodiment can 
also be represented by the magnitude of the friction circle. 
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